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Simplex Spline Bases on the Powell-Sabin 12-Split: Part II 

Tom Lyche and Georg Muntingh 


Abstract. For the space S of quintics on the Powell-Sabin 12-spht 
of a triangle, we determine the simplex splines in S and the six symmet¬ 
ric simplex spline bases that rednce to a B-spline basis on each edge, 
have a positive partition of unity, a (barycentric) Marsden identity, 
and domain points with an intnitive control net. We provide a quasi- 
interpolant with approximation order 6 and a Lagrange interpolant at 
the domain points. The latter can be used to show that each basis 
is stable in the Loo norm, which yields an bound for the distance 
between the Bezier ordinates and the values of the spline at the corre¬ 
sponding domain points. Finally, for one of these bases we provide G°, 
C^, and conditions on the control points of two splines on adjacent 
macrotriangles, and a conversion to the Hermite nodal basis. 


Analogous to the quadratic simplex spline basis from we derive quintic 
simplex spline bases on the Powell-Sabin 12-split A of a triangle A [5]. The result¬ 
ing computations are implemented in a Sage worksheet, which can be downloaded 
and tried out online in SageMathCloud [3]. We follow the notation in Part I. 

A case-by-case analysis of the possible knot multiplicities yields: 

Theorem 1. With one representative for each S3 equivalence class, these are the 
quintic simplex splines on A that reduce to a B-spline on the boundary of A; 


We hrst create a large list of potential bases for the space iSKA) of quintics 
on the 12-split. Using the macro-element from [2], we then narrow this down to a 
short list with good properties: 

Theorem 2. There are precisely six sets B — Ba, Bt, Be, Bd, Be, Bf satisfying: 

(1) B is a basis o/iS|(A) consisting of simplex splines. 

(2) B is Ss-invariant. 

(3) B reduces to a B-spline basis on the boundary. 

(4) B has a positive partition of unity and a Marsden identity, for which the 
dual polynomials have only real linear factors. 

(5) B has all its domain points inside the macro triangle A, with precisely 8 
domain points on each edge of /\. 

For instance, the basis Be = is 
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and satisfies the barycentric Marsden identity 



Factoring the dual polynomials and replacing ‘c^’ by ^Vi, one obtains 39 sets 
{Pyj-}r=iJ J = Ij • ■ ■; 39, of dual points. Taking the average of each set one arrives 
at the domain points To preserve the symmetry of A, the domain points 

are forced to form a hybrid mesh with triangles, quadrilaterals, and a hexagon in 
the center. This mesh is shown below on two adjacent macro triangles, together 
with an ordering of the domain points. 



Vl 


The collocation matrix is nonsingular, showing that is uni¬ 

solvent for iS|(A), i.e., there is a unique Lagrange interpolant at the domain points. 
Moreover, it was previously shown that there is a unique Hermite interpolant for 
the space iSKA) based on values and derivatives at the corners, midpoints, and 
quarterpoints [2]. Finally, the Marsden identity yields that 


39 5 .. 
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Pi.n + 


' Pi,rk 


is a quasi-interpolant that reproduces all polynomials up to degree 5 and has ap¬ 
proximation order 6. Moreover, using the Lagrange interpolant we show that the 
six bases are stable in the Loo norm with a condition number bounded indepen¬ 
dently of the geometry. As a consequence we obtain an bound of the distance 
between the Bezier ordinates and the values of the spline at the corresponding 
domain points. 
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As in the above figure, let A := [t’l, ^ 2 ,'^ 3 ] and A := be triangles 

sharing the edge [i>i, V 2 ]- Imposing a smooth join along [vx,V 2 \ of 


39 

I{'v) ■= X! v G A, 

2=1 


39 

/(■w) := X! v G A 

i=l 


translates into linear relations among the Bezier ordinates Ci and Ci. 


Theorem 3 . Let (/3 ij/32 , A) be the barycentric coordinates of with respeet to 
the triangle A. Then f and f meet with 
C® smoothness if and only if Ci = Ci, for i = 1 ,... ,8; 
smoothness if and only if in addition 

C9 = PlCi + P2C2 + PsCg, Cii = / 3 i( 2 c 3 — C2) + /I2C4 + fdscii, 

Cio = P1C2 + P2C3 + AciO) C12 = + / 33 C 12 ) 

and analogous conditions for Ci3,ci4, and 615; 
smoothness if and only if in addition 

CI6 = /^iCl + 2 / 3 i/ 32 C 2 + / 3 |c 3 + 2 / 3 i/ 33 Cg + 2/32/33C10 + P3C1Q, 


Ci7 = fifC2 + I3ic4 + Picn + 2l3ij32 


3c3 —C2 


■ 2/3i/33 


3cio —C 2 


2/32/33 


Ci8 = /32 2 c3±|U^ + /3|£4H^ 2A/32 ^^~^°=^^^""+°" 

I 2/3i/3^ 9 ci 2—2C5—Cll 


cio+2cii—C 3 

2 • 


and analogous conditions for Ci 9 ,C 2 o, aud C 21 . 


Whenever the domain points follow the shape of the macro triangles, we recover 
the classical Bezier conditions. All conditions are valid for the domain points as 
well, so that they also hold for the control points. Although conditions for 
smoothness can also be derived, one of these involves only (/3i,/32,/33) and the 
Bezier ordinates on one triangle, showing that this element cannot be used to 
obtain smoothness on a general triangulation. 

One can easily convert between and the Hermite nodal basis from [5]. For 
instance, the nodal function corresponding to the point evaluation at Vi is 



which, on a regular hexagon split at its barycenter, has the graph and wireframe 
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